A unified description of static and dynamic properties of Fermi liquids 
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In Landau's phenomenological Fermi- liquid theory (FLT), most physical quantities are derived from 
the knowledge of the energy variation 8E[Sn] corresponding to a change 8n of the quasi-particle (QP) 
distribution function n = {nk CT }. We show that the internal energy E[n\ (or, more precisely, the 
thermodynamic potential $[n]), expressed as a function of the QP distribution n, can be interpreted 
as an effective potential (in the sense of field theory), which is obtained from the free energy by a 
Legendre transformation. This allows to obtain explicitly <5<3> (or SE) starting from a microscopic 
Hamiltonian and to relate the Landau / function to the forward-scattering two-particle vertex 
without considering the collective modes as in the standard diagrammatic derivation of FLT. Out- 
of-equilibrium properties are obtained by extending the definition of the effective potential to space- 
and time-dependent configurations. $[n] is then a functional of the Wigner distribution function 
n = {nu.cr(r,t)}. It contains information about both the static and dynamic properties of the Fermi 
liquid. In particular, it yields the quantum Boltzmann equation satisfied by nk CT (r, t). Finally, we 
show how 5<6[<5n] can be derived (in the static case) using a finite-temperature renormalization-group 
approach. In agreement with previous results based on this technique, we find that the Landau / 
function is defined by the fixed-point value of the O-limit of the forward-scattering two-particle 
vertex. 
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PACS Numbers: 05.30.Fk, 71.10.Ay, 71.10.Ca 
I. INTRODUCTION 

Landau's original approach to Fermi-liquid theory 
(FLT) is phenomenologicallH The main assumption is 
the existence of a one-to-one correspondence between 
the low-energy elementary excitations (quasi-particles 
(QP's)) of the Fermi liquid and the elementary exci- 
tations ('particles' and 'holes') of the non-interacting 
fermion gas. Thermodynamic quantities are derived from 
the knowledge of the energy variation SE[Sn] correspond- 
ing to a change Sn of the QP distribution function n = 
{n-ka}- [k and a denote the QP momentum and spin.] 
6E[Sn] is parameterized by the Landau / function (or, 
equivalently, the Landau parameters) which describes the 
interaction between two QP's. Non-equilibrium proper- 
ties are obtained by generalizing the energy functional 
SE[Sn] to space- and time-dependent configurations, and 
by simultaneously considering the Boltzmann transport 
equation satisfied by the QP distribution. 

Much effort has been dedicated to justify Landau's 
phenomenological FLT from a microscopic Hamiltonian. 
The standard derivation consists in showing that the QP 
dynamics obtained from the microscopic two-particle ver- 
tex agpees with the conclusion of the phenomenological 
FLT.Bu Although this demonstration provides a micro- 
scopic definition of the Landau / function, it does not aim 
at calculating explicitly the energy functional 5E[5n\. 

The so-called statistical FLT provides an alternative 
approach which is in spirit much closer to the phe- 
nomenological theory of Landau. The main idea is to 
express the thermodynamic potential $[n] obtained from 
a microscopic Hamiltonian as a function of QP occupa- 
tion numbers (see Refs. 0-o and references therein). The 



equilibrium QP distribution function is determined by 
minimizing the thermodynamic potential. The excita- 
tion energies are denned as functional derivatives of the 
energy with respect to the distribution function. They 
are usually referred to as statistical QP energies to dis- 
tinguish them from the QP energies obtainod from the 
poles of the single-particle Green's function.B Although 
this formulation provides a natural bridge between mi- 
croscopic models and phenomenological FLT, few explicit 
calculations have been done in this framework. For in- 
stance, the Landau / function (which can be obtained 
as the second functional derivative of the energy with 
respect to the QP distribution function) has not been 
related to the forward-scattering twosparticle vertex as 
was done within the standard FLT.mj Besides, it is not 
clear how the collective modes of the Fermi liquid can be 
studied within this approach. 

In this paper we propose a unified description of static 
and dynamic properties of Fermi liquids. The main idea 
is to interpret the thermodynamic potential 3>[n], ex- 
pressed in terms of the QP distribution function n, as 
an effective potential (in the sense of field theory) , which 
is obtained from the free energy by a Legendre trans- 
formation. A similar definition of the thermodynamic 
potential <&[n] can be found in the statistical FLT and in 
particular in the work of Balian and De DominicisD This 
allows to obtain explicitly <5<f> (or SE) starting from a mi- 
croscopic Hamiltonian and to relate the Landau / func- 
tion to the forward-scattering two-particle vertex with- 
out considering the collective modes as in the standard 
diagrammatic derivation of FLT. Dynamic properties are 
obtained by extending the definition of the effective po- 
tential to space- and time-dependent configurations. $[n] 
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is then a functional of the Wigner distribution function 
n = {n.ko-(r, t)}. We show how we can extract from $ 
both the response functions of the Fermi liquid and the 
quantum Boltzmann equation satisfied by n-k a (Y,t). 



into account. We use a grand canonical formalism at fi- 
nite temperature T, the limit T — > being taken at the 
end of the calculations. We take h = ks = 1 throughout 
the paper. 



Outline of the paper 



II. PHENOMENOLOGICAL FLT 



In section |fij, we briefly review some basic results of the 
phenomenological FLT that are useful in subsequent sec- 
tions. The microscopic definition of the thermodynamic 
potential $[n] is given in section [II. We first consider 
the free energy in the presence of an external source field 
/ikcr that couples to the QP number operator. The QP 
distribution function n is then obtained by taking the 
functional derivative of the free energy with respect to 
the source field. By performing a Legendre transforma- 
tion, we obtain a functional $[n] of the QP distribution 
n. In field theory, $[n] is known as an effective poten- 
tial. The energy variation SE[Sn], or equivalently the 
Landau / function, is obtained from the variation <5$ [Sn] 
of the thermodynamic potential <fr[n] around its equilib- 
rium value. As a first example, we calculate exactly 
for non-interacting fermions. In the case of interacting 
fermions, we use general properties of Legendre trans- 
forms to express <5$[<5n] as a function of the (linear) re- 
sponse function Xaa' (k, k') to the external field h kcr . We 
give a microscopic definition of the Landau / function in 
terms of the susceptibility \- Using the usual assump- 
tions of FLT, we calculate x an d recover the standard 
microscopic definition of /. 

In section IV, we extend the definition of the effec- 
tive potential $[n] to space- and time-dependent config- 
urations by considering the Wigner distribution function 
ri = {nk CT (r, t)}. Restricting ourselves to small deviations 
<Jnic ff (r,t) from equilibrium, we are able to compute the 
corresponding variation of the effective potential. 

In the static and uniform limit (^(r, t) — n^), 5& re- 
duces to the Landau's functional studied in section [II. 
In out-of-equilibrium cases, by requiring (5$ to be station- 
ary with respect to variations of the Wigner distribution 
function <5nkcr(r, t), we obtain the quantum Boltzmann 
equation satisfied by Jnk<r(r, i). We also show how the 
charge and spin response functions are related to the ef- 
fective potential. 

In section [v|, we use a finite-temperature renormaliza- 
tion-group (RG) approach to reproduce the results of sec- 
tion [IL This approach provides an alternative method 
for computing the Landau's functional <5$[(5n] ancLcon- 
firms previous results based on this technique. Dill] We 
find that the / function is given by the fixed-point value 
of the fi-limit of the forward-scattering two-particle ver- 
tex (r°*). 

Unless otherwise specified, we consider a three- 
dimensional spin-i fermion gas (of volume v) with a 
spherical Fermi surface, and assume isotropy in spin 
space. Only short-range repulsive interactions are taken 



In this section, we summarize the basic results of the 
phenomenological FLT which are useful for our purpose. 



A. Landau's functional S&[5n] 

In the ground state of a Fermi liquid, the QP distri- 
bution function corresponds to a filled Fermi sea: = 
0(/cf — fc), where kp is the Fermi momentum and Q(x) 
the step function. In the phenomenological FLT, one 
postulates that a change Sn^ — - 
distribution induces the energy variatio 



£1 of the QP 



8E[5n] = ^ £ k<5n kc r 

k,<7 

+7T /^'(k.kO^nwfekv, (2.1) 

k,k',<r,o-' 

neglecting terms of order 0(6n 3 ). ek is the energy of a QP 
in the absence of other excited QP's (for simplicity, we 
consider only cases where ek does not depend on spin) . In 
the vicinity of the Fermi surface, it can be written as ek — 
v* F (k — hp) + /i, where v F is the QP Fermi velocity and 
\i the chemical potential. The term of order 8n 2 in (2.1) 
comes from the interaction between quasi-particles. For 
states near the Fermi surface, k ~ Uf and k' ~ fci?, the 
function / CTO -'(k, k') = faa'id) depends only on the angle 
6 between k and k'. The Landau parameters F t a and 
F ; a are defined by expanding / on the basis of Legendre 
polynomials: 



2iV(O)jw(0) 



J2(F?+aa'FnPi(cose), (2.2) 



where N(0) = k F /2it 2 v F is the density of states per spin 
at the Fermi level. 

In the grand-canonical ensemble at finite temperature, 
the physical quantity of interest is the thermodynamic 
potential $[ra] = E[n] - fiN[n] - /3~ 1 S{n}. N[n] = 
Sk cr rikcr ^ S ^ ne total QP number, S[n] the entropy, and 
j3 = 1/T the inverse temperature. Because of the cor- 
respondence between excitations of the non-interacting 
fermion gas and QP excitations, the cxr^'opy has the same 
expression as in a perfect Fermi gas:l 



S[n] 



][nua hin ko - + (1 - n kc r) ln(l - ri ko -)]. (2.3) 
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At equilibrium, the QP distribution function n is ob- 
tained from the condition 5^[n]/Snu.a-\n — 0: 

nk CT = n F (e k - ^i), (2.4) 
where np(x) — (e^ + l) -1 is the Fermi-Dirac factor, and 
SE\n] 



Sn ko 



(2.5) 



k',o-' 



is the QP energy corresponding to the equilibrium dis- 
tribution function n. Thus, if we expand $[n + Sn] = 
$[n] + <5$[^n] around its equilibrium value, we obtain to 
lowest order in Sn 



k.k',er,a-' r 



k,k' 



n' F (e k - /i) 



-/ CTCT <(k,k') 



5n ka 5n k ,, 



(2.6) 



where n' F (x) — — /3/4 cosh 2 (/3x/2). There is no linear 
term in Sn since is stationary at equilibrium. 



B. Deformation of the Fermi surface 



Because of the thermal factor l/n' F (ek — At) in (2.6), 
small variations of the thermodynamic potential corre- 
spond to QP excitations lying in the thermal broadening 
of the Fermi surface (|£k — /x| < T). When T — > 0, these 
excitations have vanishing energies and can be viewed as 
resulting from a displacement u CT (k) of the Fermi surface. 
Here k = k/fc is a unit vector in the direction of k. It is 
sometimes comaenient to express <5<& directly in terms of 
w CT (k). WritingEl 



Sriko- = v F u a (k)n' F (e k - u), 



(2.7) 



we obtain 



+N(0) 



47T 47T 



(2.8) 



in the limit T — > 0. We have used ek — > ek and 
n' F (x) —y —S(x) when T — > 0. fig. denotes the solid an- 
gle in the direction of k. Eq. ( plf ) was first derived by 
Pomeranchuk considering the change in energy SE result-, 
ing from a Fermi surface displacement u a (k) at T = O.tHl 
By requiring <&[u] to be minimum for u = 0, we easily 
deduce from (2.8) the stability conditions-for a 3D Fermi 
liquid: F? > -21 - 1 and F? > -21 - 1.E3 



C. Dynamic properties 

The Landau's functional <5$[<5n] allows to compute the 
thermodynamic properties of the Fermi liquid, but does 
not contain any information about the QP dynamics. To 
obtain the latter, one has to extend J-he definition of 5E 
to out-of-equilibrium configurations :Bfl 



SE\Sn c 



J2 J d 3 re k Snt(r,t) 

k,er 

+ ^ E J d 3 rf ^(k,k<) 

kjk'jfTjCr' 

x5n&(r,t)6n$ a ,(T,t). (2.9) 



n kcr( r ' *) = n k°i + ^ n k(j( r ' *) gi yes the probability of find- 
ing a QP with momentum k and spin a at point r in 
space and time t. We use the notation n cl to emphasize 
that this approach is semiclassical, since it is assumed 
that one can simultaneously specify the momentum and 
position of the QP. For low-energy excitat ions , the in- 
teraction between QP's (last term of Eq. (2.9)) can be 
assumed to be local in space. The time-dependent QP 
energy is defined by 



£k(r,f) 



SE\n 



ell 



(2.10) 



k' a' 



Eq. (2.£) is_siipplemented with the Boltzmann trans- 
port equatiorocl 



dt 



- VkO,*)- V r e k (r,t) 
+V r n k 1 CT (r,t)-Vkek(r,t)=0, 



(2.11) 



where we have used the quasi-classical equations dr/dt = 
Vk£k(r, 1) and dk/dt = — V r £k(r, t). To first order in 
Sn cl , Eq. ( 2.11 ) reduces to 



dt 



v* F k-V r 5nt(r,t) 



-5(e k -/V jw(k,k')<^k • V r H', ff ,(r,t) 

7/ fc * 



= 0. 



k',cr' 



The solution of this equation can be written as 
Snt(r, t) = v* F u a (k, q, Sl)5(e u - M )e^"- nt ) 



(2.12) 



(2.13) 



where u a (k, q, fi) is a dynamic displ acement of the Fermi 
surface. The solutions of Eq. ( 2.1 2| ) correspond to zero- 
sound (itf = itjj and spin- wave («j = — itj.) modes. 
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III. LANDAU'S FUNCTIONAL $ AS AN 
EFFECTIVE POTENTIAL 

In statistical mechanics and in the theory of phase 
transitions, it is natural to consider the Legendre trans- 
form of the free energy. In magnetism for instance, the 
Gibbs thermodynamic potential, which is a functional of 
the magnetization M(r), is obtained from the_free en- 
ergy by performing a Legendre transformation. Ej In field 
theory, such. Legendre transforms are known as effective 
potentials O In this section we show that the Landau's 
functional <I>[n] can be naturally interpreted as an effec- 
tive potential. 



A. Effective potential $ 

In order to define the effective potential we first cal- 
culate the free energy — In Z[h] in presence of a (real) 
external source field that couples to the QP occupa- 
tion number operator n^. Then we obtain the Landau's 
functional $[?i] by performing a Legendre transforma- 
tion. 

We write the partition function in the Matsubara for- 
malism as a functional integral over Grassmann variables 
ip^ (k, r) (r is an imaginary time): 



Z[h) = / V^*V^e 



(3.1) 



where S[ip*,ip] is the action when /i k(T — 0. The source 
field contributes to the action a term 



E 



hka / dr nkcr (t) 
Jo 



(3.2) 



Note that we do not specify at this point the expres- 
sion of nko-(T) as a function of the "0 field. Only for 
non-interacting fermions do we have the simple relation 
ftka-(T) = ip*(k.,T)ip (r (k.,T). The QP occupation number 
"kcr = (^ka-) is obtained from the functional derivative of 
the free energy: 



1 5 In Z[h] 



5h^ a 

Now we introduce the Legendre transformE! 
$[n] = — — In Z[h] - /i kcr n k(T , 



ft 



k,<r 



(3.3) 



(3.4) 



where /i kcr [n] is obtained by inverting ( [3.3] ) . The equation 
of state for the thermodynamic potential <5[n] reads 



8n k(7 



(3.5) 



At equilibrium, i.e., in the absence of source field, <J>[n] 
is stationary with respect to small variations of the QP 
distribution function. 



B. Non-interacting fermions 

As a first example, we consider non-interacting 
fermions. Using hk, a (r) — ^/£(k, r)^ CT (k, r), we obtain: 

Z[h] = J Vijj*ViPe^ k ^ rAK " ){lu '~' ik ~ hk " +t * ),f ' Ak ' uj) 

=n( i+e ~ ,3(ek+ ' ik ^ i) )' ( 3 - 6 ) 



k.cr 



where e k is the energy of a fermion with momentum 
k. We have introduced the Fourier transformed field 
tp^\\s.,u>) where lu = irT(2m + 1) (m integer) is a 



fermionic Matsubara frequency. Eq. (3.3) yields 

Hka = n F (e k + h ka - fx). (3.7) 
Inverting ( |3.7| ) and using ( |3.6[ ), we eventually obtain 

* N = Vl^k - A 1 )"-kcr + -3 VVk<x hinka 

(3.8) 



k,cr k.cr 

+ (1 - n-krr) ln(l - n kcr )]. 



Eq. (3.5) is the expected result for non-interacting 
fermions. 



C. Interacting fermions 

For interacting fermions, it is not possible to calculate 
exactly the thermodynamic potential <&[n]. However, we 
do not require the whole knowledge of $[n], but only 
its variation <5$[5n] when the QP distribution function n 
varies from its equilibrium value n = n\h=o by an amount 
Sn. [For T — > 0, n kcr = 0(kp — k) corresponds to the 
ground-state distribution function.] This turns out to be 
a much easier task than obtaining $[n]. 

Expanding $[n + Sn] to second order in Sn, we obtain 



= 5 E 



«J(2)$[ 



2 , h* 1 , SfDarSnVai 
k.k ,rr,c7 



<5n ko .(5n k /, 



(3.9) 



There is no lin ear term since $ [n] is stationary at equilib- 
rium [Eq. (3.5)]. To proceed further, we use the following 
relation between the functional derivatives of the free en- 
ergy r=j/3 _1 In Z[h] and those of the Legendre transform 



E 



Z?- 1 ^ 2 ) \nZ[h] 



£(2)$[ 



h=0 



Snk 3 a 3 Snk 2Cr2 



Introducing the matrix 



&<T\ ,(72 ^ki ,k2 ' 

(3.10) 
(3.11) 



h=0 
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we thus obtain 

= 5 E X^(k,k')5n kff <5n k v. (3.12) 

k, k ',<7,<T' 

Note that Xctct' (k, k') is nothing but the (linear) response 
f unc tion t o th e external source field /i k(T . Comparing 
@ and (PI ) (and assuming that the system is a Fermi 
liquid), we obtain a microscopic definition of the Landau 
f function: 



1 



faa' (k, k ) , 



" W,Vk ' k ' ' X^(k,k'). (3.13) 



From now on, we ignore the difference between e k and 
e k since we are ultimately interested in the limit T — *■ 0. 
Thus, the calculation of the Landau / function reduces 
to the calculation of the susceptibility x- I n the next sec- 
tion, we show how \ can be calculated using the standard 
assumptions of FLT. 



D. Calculation of the Landau / function 

In this section, we first express the Landau / function 
in terms of the QP properties, and then relate it to the 
bare fermion properties. 



1. f in terms of QP properties 

Since the external field /i k(T couples to the QP's, we 
have to distinguish between the QP field (that we de- 
note by ip*)) and the field corresponding to the 
bare fermions. With these notations, we have n k(T (T) = 
tp* (k, T)Vy( k, t) , where t he o perat or n ka - (r) has been in- 
troduced in ( |3.2| ). From ( |3.2|) and ( 3.1l| ), we obtain 



(3.14) 



where x is the two-particle Green's function 

Xaa>(k,k';q) = (tiCk-)Mk+)^CK)^'C k '-))c 
= (^Ck-)Mk+)^Ck+)^4k^ 

-^,o(^(fc)^(fc))(^(fc')^(fc'))- 

(3.15) 

The average values are evaluated in the absence of source 
field (h = 0), and the meaning of (• ■ -) c is defined by the 
second line of ( 3.15 ). We use the notation k — (k,iuj), 
q = (q, ifl), and k± = k±q/2. lo and Q denote fermionic 
and bosonic Matsubara frequencies, respectively. The 
Landau / functio n is related to the limit q — > of 
Xaa* (k, k'; q) [Eq. ( 3.13j )] . In a Fermi liquid, the forward- 
scattering limit g — > is ill-defined and one should distin- 
guish between the q- and fJ-limits, which correspond to 



rZ/q — > and g/fi — > 0, respectively. Since a static and 
uniform external field cannot create quasi-particle-quasi- 
hole excitations (and therefore modify the ground-state 
QP distribution function) in the fl- limit ,E£l / is obtained 
from the g-limit of the function x- It turns out to be 
more convenient to keep a finite q at intermediate stages, 
the limit q — > (with 0,/q — » 0) being taken only at the 
end of the calculations. Besides, the function x (with q 
finite) will also have to be considered in the analysis of 
the dynamic properties of Fermi liquids (section [TV]) . 
X satisfies the Dyson equation 

Xaa'(k,k';q) = Xav>(k,k'-,q) 

E Aikk-A) 

k 1 ,k 2 ,o\ ,ct 2 

xf '^(h, h\ $)w(*a, k';q), (3.16) 

where 

f %,{k, k'; q) ee T™,^ a (k + ;k'_~k' + ,k-) (3.17) 

is the irreducible QP vertex in the Landau channel.0 
The non-interacting part 

X ( °l(k,k';q) = -8 a ,„'5 % #G(k+)G(k-) (3.18) 
is easily expressed in terms of the QP propagators 



G(k±) = -(Mk±WAk±)) 

...n q 

ilu ± i— - e k =F v k • - 



(3.19) 



where e k is the QP energy and v k = V k e k the QP 
group velocity. Approximating e k by v* F {k — kp) + jjL, 
we have v k = v F k (with k = k/fc). The quantity x' ' 
becomes singular in the forward-scattering limit q — » 
since the poles of the two QP propagators coalesce. In 
FLT, one assumes that x^ is the only singular quantity 
in this limit This implies that the irreducible vertex 
r irr , which does not contain X , is a non-singular quan- 
tity that has a well-defined limit when q — > 0. We there- 
fore set q = in f lrr . Furthermore, since the singularity 
of x^ is due to QP states in the vicinity of the Fermi sur- 
face, we can ignore the k and uj dependence of r irr (i.e., 
set k = k' = kp and u) = lo' — 0), which then becomes 
a function of k and k'. [The variables \k — kp\,\k' — kp\ 
and w,u>' are irrelevant in the RG sense (see section^).] 
This allows to perform the frequency sums in ( 3.16| ) and 
obtain a Dyson equation for x : 

W (k,k';g)= X ^,(k,k';g) 



This can be rewritten as 



\ J2 A(*M;q) 

kl ,k2 ,<Tl ,<T2 

(k b k 2 ) X „ 2 y(k 2 ,k'; ? ~). (3.20) 
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X^(k,k';g) = ^(k,]^) + ^-0(k,k'). (3.21) 



In Eq. X _1 (9)j X ( ° Kl (g) and r irr should be under- 

stood as ma trices in sp in (a ) and momentum (k) space. 
Comparing ( 3.13 ) and ( 3.21 ), and noting that 



xi J'( k > k '; ( ?) = 5<T,<T'S±,vri F (eu.- (i) 



Vk • q 



iQ — Vk • q 
Sa^'Swn'ptek - n) (3.22) 



in the g-limit, we conclude that 

/ CTCT ,(k,k')=O(k ) k0. 



(3.23) 



We can proceed one step further by relating r irr to the 
total QP vertex F, using the equation 

f aa , (k, ~k'; q) = r£, (k, k') - i- ^ r£ t (k, kx) 

fcl ,&2 :°"2 

xx^U (fci, fc 2 ; ?)f ^(fc, fc'; g). (3.24) 

Since F, as F lrr , is independent of frequencies, we can 
perform the sum over cjj and Cl>2 in Q3.24j ), which yields 

t aa , (k, k' ; $) = r£, (k, k') - i ^ f ^ (k, k x ) 

kl ,k2 ,<Tl ,(72 

xx^^kx,^; g)f CT3 ^(k 2 ,k'; g). (3.25) 



X^ -* vanishing in the fi-limit [Eq. ( 3.22| )], we obtain 



O(k,k0=^im[f CTCT ,(k,k';g) 

_ fSf2 



9=0. 



r"„'(k,k'). 



(3.26) 



Eqs. ( 3.25 , 3.26 ) relate the Landau / function to the Sl- 
limit of the forward-scattering QP vertex. 



2. f in terms of bare fermion properties 

The last step of our derivation is to relate r irr to the 
vertex of the bare fermions. Consider first the single- 
particle propagator in a Fermi liquid. It can be written 



G(k) = -(w(fe)C(fc)) 



ILJ — Ck + fJ. 



+ Gi nc (fc), 



(3.27) 



where Gi nc (£;) is the incoherent part, and z the QP renor- 
malization factor. The QP field can be obtained (at 
least formally) by filtering out the incoherent part and 
rescaling the field in order to ensure that the total spec- 
tral weight is equal to one (see Nozieres' bookta for a 
detailed discussion). [As discussed in section [yj the RG 
approach provides an explicit realization of this proce- 
dure.] G and G are then related by 



G{k) = zG{k) + G iac (k). (3.28) 

Consider now the two-particle Green's function 

G" a ,(k,k';q) = mk-)Mh)^(k' + )^'(k'-)} c . 

(3.29) 

G u can be written as 
Gt,(k,k';q) = G I W\k,~k'- ) q) 

£ G^ihh-A) 

k 1 ,k 2 ,cri,a2 



xr„ 



(h,h;q)G^J(k 2 ,k. f ;q), (3.30) 



where 



r aa' {k, k ; q) — r<T<T',(T',(T (fc+ 7 k_ , fc^_ , k— ) (3.31) 
is the total vertex for the bare fermions, and 

G 1 S\~k,k';q) = -5 a , a ,5 Wk ,G{k + )G(k-) 

= -8^,8 k y[z 2 G{k + )G{k.) + V {k)). 

(3.32) 

We have separated the coherent part of the particle-hole 
pair propagation from the incoherent part (p(k). Note 
that ip(k) is non-singular in the lim it q — > an d is t hcrc- 
fore evaluated at q = 0. Inserting ( 3.32| ) into ( 3.30 ) and 
retaining only the coherent part of the particle- hole prop- 
agation, we obtain an equation for the Green's function 



X introduced in section III D 2 
Xaa>(k,k';q) = x°l> (k, k'; q) 



kl.k 2 ,(T1 ,<T2 

xz 2 v airT2 (h, k; §)x£U(fe, ?)■ (3-33) 



We then deduce from J3.16| ) and ( [3.33] ) 

f aa i (k, k'; q) = z 2 T aa , (k, k';q). 



(3.34) 



Eq. (3.34) relat es the QP vert ex F to the bare fermion 
vertex r. From we obtain the well-known 

expression of the Landau / function in a Fermi liquid: 



/ ffff ,(k,k') = z 2 r^(k,k') 



(3.35) 



IV. DYNAMIC PROPERTIES 

The Landau's functional $[n] does not contain any in- 
formation about the dynamic properties of the Fermi liq- 
uid. In this section, we extend the definition of the ef- 
fective potential <E> to space- and time-dependent config- 
urations by considering the Wigner distribution function 
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(WDF) n = {n ka (r,t)}. This allows to obtain the QP 
dynamics and the collective modes of the Fermi liquid 
without introducing the semiclassical distribution func- 
tion H kCT (r, i) and the corresponding Boltzmann trans- 
port equation as in the phenomenological FLT (see sec- 
tion flQ . 



A. Generalized effective potential $ 

The best quantum analog to the semiclassical distribu- 
tion function njj^r, t) is the WDFEj The latter is not a 
true distribution function since it is not positive definite. 
However, as far as its moments are .concerned, it behaves 
similarly to a distribution function.Eao For our purpose, 
we define the WDF in Matsubara time as 



n k(T (r, t) 



d 3 r'e~ 



X (K{ T -T T )H T + 1>' 



(4.1) 



and its Fourier transform as 



\- f dT I d 3 re- 4q r+ ^ T n kff (r,r) 
™ Jo J 



~£(C(fc-)Vv(fc+)>. 



(4.2) 



Note that these definitions involve the QP field 

We are now in a position to define an effective potential 
which is a functional of the WDF n = {n k(T (<?)}. 
We shall proceed along the same lines as in section 



[II. We consider the system in presence of a source 
field ft. ko -(<z) = h^ a {—q) that couples to the QP operator 
nk a (q) = P^ 1 Y],,, i >g (fc_)Vv(fc+). We write the partition 
function as in (pTj), with 



k,tr,g 



(4.3) 



The WDF is obtained by taking the functional derivative 
of the free energy with respect to the source field: 



n-ka{q) 



1 SlnZ[h] 
(n^aiq)) = - -r- 



Sh ka (-q) 
The effective potential is then defined as 

$[n] = ~lnZ[h] - h ka (-q)n ka (q), 

and satisfies the 'equation of state' 
5$\n] 



Sn kcr (q) 



rr = -h ka (-q). 



(4.4) 



(4.5) 



(4.6) 



We will show below that in the absence of source field 
(h = 0), the stationarity condition of the effective po- 
tential [Eq. ( |4.6| )] yields the equation of motion for the 
WDF. 



Even for non- interacting fermions, $[n] cannot be cal- 
culated exactly. We shall therefore consider only small 
fluctuations <5n ko .(<7) around the equilibrium state: 



n-kaiq) = Sq^n^a + Sn ka -(q). 
Eq. ( 3.1 0[ ) can be easily generalized into 



(4.7) 



5n k<T (-q)Sn k i a >(q') 



= hs'X-^y-,q), (4.8) 



where 



Xaa'Q*-, k';q) 



prH& hxZ[h] 



Sh ka (-q)5h k/l7l (q) 



(4.9) 



h=0 



is the susceptibility introduced in section III 
[Eqs. (|T4|j3T|)]. The Kronecker symbol 5 M , in ( pf 
results from translational invariance. To second order in 
5n, we thus have 



5<f>[6r 



X a l> (k, k'; ^Sn^i-^Snva, (q). 



(4.10) 



Using Eqs. ( |3 . 2 1| , |3 . 22|j3 . 23| ) , we can write the effective 
potential in terms of the Landau / function: 



1 \p f 5 a ^,5 k . k , iQ, - v k ■ q 

2 . . ^ , X n' F (eu - (J,) v k • q 



+ -U a , (k, k')\6n ka {-q)Sn k , a , (q). (4.11) 
v > 



Eq. (4.11) generalizes the Landau's functional to space- 
and time-dependent configurati ons . As in the static and 
uniform case studied in section III, <5$[<5n] is essentially 
parameterized by the Landau / function. We shall show 
in the next section that <5$ contains all the information 
about the QP dynamics. Moreover, if we consider the 
q-limit, the WDF {n ka .(r, t)} reduces to the QP distri- 
bution function {n kcr }, and <&\n\ to the thermodynamic 
potential introduced in section III . Thus the effective po- 



tential $ describes static and dynamic properties of the 
Fermi liquid in a unified framework. This result should 
be contrasted with the phenomenological FLT which re- 
quires slightly different approaches to deal with static and 
dynamic properties (see section ||), and strongly relies on 
semiclassical arguments for the latter. 



B. Quantum Boltzmann equation 

In the absence of source field (h = 0), the stationarity 
condition of the effective potential [Eq. ([1.6])] yields 

(Vk • q - ifl)Sn ka (q) - v k • qri^(e k - n) 

x - Y jw (k, k')6n k , a , (q) = 0. (4.12) 

7/ Z J 
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Eq. (4.12) is the quantum Boltzmann equation for the 
WDF. [Real-time quantities are obtained by the usual an- 
alytic continuation iQ — > fi + i0 + .] Note that in the sim- 
ple case we are considering here (no external electric field, 
small fluctuations around the equilibri um s tate, etc), it 
is identical to the Boltzmann equation ( 2.12 ) satisfied by 
the semiclassical d istrib ution function n^ a (r, t) ll3'E2l 
The solution of ( 4 . 1 2j ) can be written as 



<5n k <r(<?) = VpUaQs., q)n' F (e k - (J,), 



(4.13) 



where ?v(k, q) is naturally interpreted as the dynamic 
Fermi surface displacement for spin a fermions. It satis- 
fies the equation 



-v k • qN(Q) 



E 



dCL f 



4tt 



(v k • q - ifl)u a (k, q) 
i/^(k,k'K'(k',?)=0. (4.14) 



Here we h ave used n' F (e-^ — fi) = — <5(ek — /i) when T — ► 0. 
Eq. ( 4.14 ) yields the zero-sound (for W[ = u{) and spin- 
wave (for it | = — itjj modes of the Fermi liquidBH 

It is also possible to express the effective potential di- 
rectly in terms of the dynamic Fermi surface displace- 
ments: 



5<f>[u] 



vv. 



2 N(0) 



E / 



4tt 



~\ |2 



Vk • q 



x|u CT (k, q)\ 2 + N{0) 
ff (k,- g >Uk',g)} 



47T 47T 



/^'(k,k') 



xu. 



(4.15) 



Here we h ave u sed it CT (k, —q) — u* a (k, q). The equations 
of motion ( 4.14 ) are then directly obtained from the sta- 
tionarity conditio n <5< l>/fa rr (k, — q) — 0. Note also that 
( 4.15 ) reduces to (|2~8[) in the q-limit. 



C. Response functions 

Since the effective potential is essentially determined 
by the susceptibility Xaa' (k, k'; q), it also yields the re- 
sponse functions of the Fermi liquid. For instance, the 
charge and spin response functions are given by 



h(q) = - E X<r<r'(k.k';g), 

k,k', 

Xs P (<?) = - E gg 'Xcjc' (k, k' ; q) , 



(4.16) 



k.k' 



where x is related to the effective potential by Eq. (4.10). 
Note that when calculating response functions, the dif- 
ference between particles and QP's can be ignored. This 
property is a consequence of the Ward identities that re- 
sult from particle-number conservation (see Ref. |2l| for a 
detailed discussion). 



Eqs. (4.1C.4.11) determine x 1 - If, for simplicity, we 
consider the case where only the Landau parameters Fq 
and Fq are non-zero, the matrix x~ X can easily .ie in- 
verted, and we obtain the well-known expressions^ 



Xch (?) = 2JV(0) 
Xsp (q) = 2N(0) 



£l ao (iQ/v F q) 
l + F^Qoiin/vpq)' 
Q,oa(iCl/v F q) 

l + F^noQiin/vpq) 1 



where 



f 1 du U X 35 + 1 

\Iqq\x) = I — =l--ln 



2 u — x 



1 



(4.17) 



(4.18) 



V. RENORMALIZATION-GROUP APPROACH 

Recently, RG techniques based on low-energy fermion 
effective actions have been applied to interacting 
fermions in dimension d > 2 by many authors (see 
Refs. I-HIH-H and references therein). The finite- 
temperature RG approachEZI was first applied tOrFermi 
liquids in dimension d > 2 by Chitov and SenechalE Con- 
trary to other works on the subject, it revealed that the 
effective interaction function in the Landau or zero-sound 
channel (particle-hole pairs at small total momentum and 
energy) does not stay marginal under the RG transfor- 
mation, since its /3-function is not identically zero. From 
the RG eaualions, the standard FLT results have been 
recovered.Qilj It has been pointed out that the bare in- 
teraction function of the low-energy fermion effective a< 
tion cannot be identified with the Landau / function 
The latter, along with other observable parameters of a 
Fermi liquid, is obtained at the fixed-point of the RG 
equations, i.e.^w-M 11 all degrees of freedom have been 
integrated outBlld The finite temperature RG approach 
has also given results that cannot be obtained within 
the standard derivation of FLT. Chitov and Senechal, 
taking into account the interferences between the zero- 
sound (ZS) and exchange (ZS') channels, have obtained 
RG equations that satisfy the Pauli principle contrary to 
the standard microscopic derivation of FLT.B Perform- 
ing a two-loop order calculation, the present author has 
obtained a non triviaL-expression of the wave-function 
renormalization .factor JI9 which was also, obtained from 
2D bosonizationE£l and Ward Identities .Eil 

In this section, we s how how the response function x 
introduced in section HI, and therefore the Landau / 
function, can be calculated using a a finite-temperature 
(Kadanoff- Wilson) RG approach. We follow the proce- 
dure used in Ref. [lO] to obtain the compressibility of a 
Fermi liquid. For simplicity we consider a 2D fermion 
gas. 

We write the partition function Z[h] as a functional 
integral over Grassmann variables [Eq. (3.1)], where, as- 
suming that the high-energy degrees of freedom have 
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been integrated out (in a functional sense), the action 
describes fermionic degrees of freedom with \k — < 
Ao -C kp. We write this low-energy effective action as 
S A ° + S A °, where 

sA ° = ~ E & (*) ( iuJ - £k + (*) 

+1^7; E E r&*,™&.k.*3.**) 



fel...fe 4 CT 4 
X^* 4 (fc 4 )C 3 (^3)^2 (fe)^! (fcl) 



(5.1) 



It will be shown below that the RG procedure ensures 
that the action is always expressed as a function of the 
QP field <ip(*>. In (5.1), the wave-vectors k satisfy \k — 
kp\ < A . The single-particle excitations are linearized 
around the Fermi surface: ek — fJ- = vp{k — kp) where vf 
is the bare Fermi velocity. The bare (antisymmetrized) 
two-particle vertex r A ° is assumed to be a non-singular 
function of its arguments. The summation over wave 
vectors is defined by 



1 



d 2 k 



kF 



fcj?+A 



dk 



kp-An 



2tt 



dO 



^ ^, (5-2) 



2tt : 



ignoring irrelevant terms at tr ee-level. 

the Landau's functional 



As shown in section [IB 



<5<I>[<5n] can also be expressed in terms of the Fermi sur- 
face displacements u CT (k) = u CT (0), or equivalently the 
'density' variations 8p a {9) = (kp /2ir)u a (9). [Here and 
in the following, we__denote by the direction of a given 
momentum k, i.e., k = (cos 6*, sin0).] For purely techni- 
cal convenience, we therefore consider an external field 
h a {9) that couples directly to the QP density operator 



, r) and write the action S h as 



(5.3) 



where 

r0 



r k F +A 



[ dTp a (9,T)=v F N(0)J2 I 
Jo - Jk 

x lim \i>* a {k + q)^(k) 

q— >0 L 



k F -A 



0=0 



(5.4) 



We take the g-limit as discussed in section [II. N(0) = 
kp/2wvp is the bare 2D density of states per spin. 

It is straightforward to show that the microscopic def- 
inition of the Landau / function [Eqs. (3.11,3.13)] be- 
comes: 



N(0) 



[6^,2*6(9) + N(0)U,(9)} = xl'HO), 



9') = 4tt 



[3- l 5^ hxZ[h] 



8K{9)5h a ,{9>) 



(5.5) 



The meaning of the notation \* is discussed below. 

The constraint to have all momenta in the shell \k — 
kp\ < Ao restricts the allowed scatterings to diffusion of 
particle-hole, or particle-particle, pairs with small total 
momentum (q < Ao).E3 Consequently, only two vertex 
functions have to be considered: the forward-scattering 
vertex function and the BCS vertex function. In the ab- 
sence of BCS instability, we can neglect the latter. We de- 
note by r^°(fci, k 2 ; q) the forward-scattering vertex. [We 
use the notation T A ° = r^j 1 ^ 4 .] q — (q,ifl) where 
q is the momentum transfer (q < Ao) and the bosonic 
Matsubara frequency fl the energy transfer, (fci, <ti) and 
(&2> 02 ) correspond to the two incoming particles, and 
(ki + q, (T4) and (k 2 — q,cr 3 ) to the two outgoing particles. 
Neglecting the irrelevant dependence on ki,2 and 071,2, 
we consider only the coupling function T A ° (ki , £2 ; q) = 
T^(9 1 ,9 2 ;q). 

The Kadanoff- Wilson RG procedure consists in suc- 
cessive partial integrations of the fermion field degrees of 
freedom in the infinitesimal momentum shell Aoe~ dt < 
\k — kp\ < Ao where dt is the RG generator and 
A(t) = Aoe _t the effective momentum cut-off at step 
t . Each partial integration is followed by a rescaling 
of radial momenta, frequencies and fields (i.e., u>' = slu, 
k' — kp = s(k — kp) and = ip with s = e dt ) in or- 
der to let the quadratic part of the action (5.1) invariant 
and to restore the initial value of the cut-off. The partial 
integration modifies the parameters of the action which 
become functions of the flow parameter t. In the fol- 
lowing, we note T the running (i.e., cut-off dependent) 
vertex (we do not write its i-dependence explicitly). For 
the purpose of our calculation, it is sufficient to consider 
the q- and il- limits of the forward-scattering vertex: 



Tl l {e x -e 2 )= Mm 1^(01, 02 5 9) 

9^0 L O=0. 



r 



h) = lim 



T ai {6 u e 2 \q) 



9 =0. 



(5.6) 



The RG process also generates corrections to the 
source field h a {9) along with higher-order terms in h. 
AL-step i, the source term in the action can be written 
asc3 (ignoring terms of order h 3 ) 



d9 dO' 
~2~k^v; 



K{9)(3 X oo'{9-9')h <7 ,{9 1 ). (5.7) 



h=0 



We do not write explicitly the dependence of and x 
on the flow parameter t. 

Note that z {h ^ and X are not physical observables, since 
they do not result from the integration over all degrees 
of freedom. Only their fixed-point values z^* and x*> 
obtained when A(t) = 0, are physical (observable) quan- 
tities. The Landau / function is therefore related to x* ■ 
[Note that x* corresponds to what we denoted by x m 
the preceding sections.] 
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A. one-loop order 

The integration of high-energy degrees of freedom 
( \k — kp\ > Ao) leading to the low-energy effective action 
( |5.l| ) will in general generate a wave-function renormal- 
ization factor z\ < 1 and an external field renormaliza- 
tion z^ h '\\ . We ign ore th ese complications which will be 
discussed in section VB . Thus, the initial conditions of 
the RG equations, besides r|A = T A °, are 

W(f)k=0. (5.8) 



The latter equation follows from the condition T <C Ao 
and the fact that the external field h creates excitations 
only in the vicinity of the Fermi surface. 

The external field renormalization at one-loop order is 
given by (see Fig. 2a in Ref. |l0| ) 

jnii 



dz 



(h) , 



E 



2tt ™" 



; £ / dk"^[G(k")f, (5.9) 



where G(k) — {iuo — vpik — hp)) 1 is the QP propagator. 

J indicates that the integration is restricted to the de- 
grees of freedom that are in the infinitesimal momentum 
shell to be integrated out. Using 



dk"[G(k")f = - 



we obtain 
dz™,(l) N(0)/3 R 



%-r- dt, (5.10) 

v p cosh pa 



dt 



cosh /3n 



E&(0i^^,^(0- (5-ii) 



We have introduced the dimensionless inverse temper- 
ature (3r = Vp/3A(t)/2, and expanded in circular har- 



monics the quantities appearing in (5.9), z^,(0) 



Yli z aJ'(l) e%W ■■■ Eq. ( |5.1l| ) is solved by introducing 

(o = 4?(o +4{ } (o, 



z {h \l) 



Ah) 

'T1 



(0 



Ah) 

Al 



(0, 



(5.12) 



and the spin symmetric (A q ) and antisymmetric (B q ) 
parts of the two-particle vertex defined by 

2JV(0)r« ( (Z) = AfS^Jw, +Bfr fflCT4 ■ r CT2(T3 , (5.13) 



where r denotes the Pauli matrices. Eq. (5.11) decouples 
into two independent equations: 



dlnz ( +\l) 
dt 

d\nz^\l) 
Jt 



Br 



cosh (3r 

Pr 
cosh 2 (3r 



-Ah 



Bf 



(5.14) 



Eqs. ( |5.14 ) have to be supplemented with the one- loop 
RG equations for the vertex functions A q and B q . As 
shown in Ref. 0, the latter can be written as 



dAl _ 


Pr 


Af, 


dAf 


dt 


cosh 2 (3r 




dt 


dBl 


Pr 


Bf- 


dB? 
1 


dt 


cosh 2 (3r 


dt 



(5.15) 



The first terms of the rhs of Eqs. ( 5.1 5| ) are the con- 
tribution of the ZS graph to the renormalization of T q . 
The contribution of the ZS' and BCS graphs is taken 
into account via the second terms on the rhs of ( 5.15 ) 
(see Fig. 1 in Ref. |l0|). Because the thermal factor 
Pr/ cosh 2 Pr is a strongly peaked function of A(t) near 
A(<) = whcn-JZj-^-i' 0, these equations have the approxi- 
mate solution: 



AKr) = 



A 



BKr) 



\ + {l-T)Af 

1 + (1 - T)Bf 



(5.16) 



for A(t) < T/v F . Af = Ap| A(t)=0 and Bf = 
-BpU(t)=o are the fixed-point values of Af and Bp. We 
have introduced the parameter r = tanh/?R and used 
tanh(/3w jF A /2) ~ 1 for T < A /v F . Since Eqs. f[l|) 
hold beyond one-loop order, we postpone their detailed 
derivation to section VB. Eqs. ( |5.14 ) show that the RG 
flow of z^ becomes significant only at small energy when 
A(t) < T/v F - This allows to insert Q5.16] ) into ( pU4j ) and 
obtain 



Ah) 



(0 



1 + (1 - t)A?* ' 
1 

r)SP* ' 



l + (l-r) 



(5.17) 



using the initial conditions z + (0U o = z — (0U o = 1- 

The renormalization of x is given by (see Fig. 2b in 
Ref. 0) 

a" J 77 

x^E/ dk"[G{k")f. (5.18) 



Expanding ( |5.18 ) in circular harmonics and using (5.10) 
we obtain 



rfwffl = JV(0) 
dr v 



E^(o4^(o. 



(5.19) 



To derive ( |5.19| ) from fl5.18| ), we have used the fact that 
z^,(ff) is an even function of 9. Thus, we obtain 
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dr 2v 



(It 



2v 



(5.20) 



The integration of ( 5.20 ) using ( 5.17 ) and the initial con- 
dition Xo-o-'COUo — yields the fixed-point matrix 



xfr(0 



xu(0 



jv(Q) 

2z/ 

^(0) 
2^ 



1 



1 + Ap* 1 + Bp* 



1 + A 



12* 



1 + Bp* 



(5.21) 



The inverse matrix is 



x*f x (0 



2AT(0) 
v 

2N(0) 



(2 + Ap* +Bp*), 



(A"* - Bp*) 



(5.22) 



From t he m icroscopic definition of the Landau / func- 
tion [Eq. fl5.5p 1, we deduce 

jw(0) = r?* vv (0), (5-23) 

a result fist obtained in Ref. [q. 

B. Beyond one-loop 

In this section, we show that higher-order loop contri- 
butions do not change Eq. ( [5.25 ) , except for a dependence 
of / on the wave-function renormalization factor. The 
RG equations can be solved exactly if we make the three 
following assumptions (see also Ref. [k], section 5): we 
assume the existence of well-defined QP's in the vicinity 
of the Fermi surface (i). Except the one-loop ZS graph, 
all the graphs are regular in the limit q — > (i.e., give 
the same contribution in the q- and f2-limits) (ii), and 
give a smooth contribution to the RG flow of various 
physical quantities (iii). [Because of the thermal factor 
/3/4cosh 2 fin, which is a strongly peaked function of A(t) 
around A(t) = for T — > 0, the ZS graph yields a singu- 
lar contribution (with respect to A(i)) to the RG flow of 
T q .] Assumption (ii) can be explicitly verified at one- loop 
order. Assumption (iii) is only 'approximate' in a sense 
that is further discussed below. 

Let us first consider the RG equations for the two- 
particle vertex. They can be written asa 



dt 



dt 



dt 



zs 

PR 



dA? 
dt ' 



ZS 



cosh {13 r) 



A\\ 



(5.24) 
(5.25) 



and a similar equation for B q . We have written explicitly 
the contribution of the one-loop ZS graph which distin- 
guishes between the q- and f2-limits. The second term on 



the rhs of (5.24) includes the contribution of other one- 
loop graphs (ZS' and BCS channels) as well as higher- 
order loop corrections. [Here we use assumption (ii) and 
the fact that the ZS graph vanishes in the 0- limit. 



It is tempting to neglect the term dAf/dt in (5.24) 
(see for instance Ref. |2^). Because of phase-space re- 
strictions (which result from momentum conservation at 
the interaction vertices), the corresponding diagrams are 
suppressed by the small parameter A^/kF if one consid- 
ers only low-energy states \k — kp\ < A (assuming that 
high-energy states (|fe — kp\ > A ) have already been 
integrated out). This property results from a frustra- 
tion of-the interfe rence s between channels in dimension 
d > 2.E3 Solving Q5.24 ) without the last term of the dig 
is equivalent to an RPA calculation in the ZS channeliM 
This RPA calculation can also be done using standard 
diagrammatic theory. Within this approximation, the 
Landau / function is naturally identified with the bare 
vertex function r A °|g =0 . [We do not need to distinguish 
between the q- and the f2- limits, since r A ° is a regular 
function of its arguments.] A major drawback of this ap- 
proximation is that the momentum scale Ao, which is not 
a physical scale in the problem, enters the definition of 
the Landau parameters in an essential way. [For a further 
discussion of this RPA approximation, see Ref. [|] 

In Refs. ||,[l(| the present author an d G. Chit ov have 
proposed another approach to solve (5.24 5.25j ), which 
bears some similarities with Landau's solution of the 
Bethe-Salpeter equation for the two-particle vertex. Be- 
sides the fact that it is not suppressed by the small pa- 
rameter Ag/kp, the ZS graph presents another interest- 
ing feature. Its contribution to the RG flow of A q and B q 
becomes singular (with respect to A(t)) at low tempera- 
ture since /3r/4 cosh 2 (3r is exponentially suppressed for 
A(i) > T. To see how this allows to solve Eqs. Q5.24| , |5.25p , 
we integrate the latter (using A q \\ = A |a ) : 



A g (t) =A?(t) - I 
Iterating fl5.26| ), we obtain 



^—AUt 
o cosh pr 



'\2 



A?(t)=AP(i) 



dt'- 



, A\\t'Y + 
o cosh /3 R 



(5.26) 



(5.27) 



According to assumption (iii), Af(t') is a smooth func- 
tion of the cut-off A(f ) when A(t') < T/v F . This im- 
plies A2{t')\x(t>)<T/v F — Ap* at low temperature. Since, 
on the other hand, /3r/4cos1i 2 /3r is strongly peak ed fo r 
A(f') < T/v F , we can replace Ap(t') in the rhs of Q5.27j ) 
by its fixed-point value Af*. The RG equation for A q 
then becomes (for A(t) <T/vf) 



A q (t) = A 



dt'- 

o cosh z [3 j 



-AHt'f 



(5.28) 



Eq. ( [5.28 ) corresponds to a decoupling of the ZS channel 
from the other channels: when A(t) > T/vf, the ZS 
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graph does not contribute to the RG flow, and both A q 
and A evolve smoothly when the cut-off decreases (with 
A q ~ A n ); when A(t) varies from ~ T/vp to 0, only 
the ZS graph gives a significant contribution (assumption 
(iii)) and drives A q towards its fixed-point value A q * (but 
does not contrib ute to the renormalization of A ). The 
solution of Eq . ( 5~28| ) and the analog equation for B q 
is given by ( 5.16 ). We now obtain an RPA-like relation 
between two fixed-point quantities: Af* = Af*/(1+Ay*). 
Both quantities are physical in the sense that they are 
obtained by integrating all the degrees of freedom. 

One should combine the RG equations of A q and B q 
with RG equations of the Fermi velocity vp(t) and the 
wave-function renormalization factor z(t). Since the ZS 
loop does not appear in self-energy corrections, Vp(t) and 
z(t) are smooth functi ons o f the cut-o ff A(t). Thus, the 
argument leading to ( 5.28 ) and ( 5.16| ) also holds when 
self-ener gy c orrections are taken into account. When 
solving ( 5.28 ), one should replace vp(t) and z(t) by their 



fixed-point values v* F and 

If assumption (iii) were to hold exactly, Eqs. ( 5.1(i| ) 
would be exact for A(t) < T/vp and T — > 0. It has been 
pointed out in Ref. [?] that the ZS' graph also gives a sin- 
gular contribution to the flow of T q {6 — 9') when the two 
incoming particles have parallel momenta (9 — 9' — ► 0). 
In this case, the ZS' loop becomes identical to the ZS loop 
(see Fig. 1 of Ref. [l(]). Thus, assumption (iii) is not exact 
and Eqs. ( 5.16| ) hold only approximately. The singular- 
ity in the ZS' channel is however restricted to very small 
angles \9 — 6'\ < T/vpkp. Only for those angles do the 
ZS and ZS' channels interfere when A(t) < T/vp. Conse- 
quently, only the components T q (l) with I > vpkp/T are 
affected by the small-angle singularity in the ZS' chan- 
nel. For most physical quantiti es (s pecific heat, effective 
mass, compressibility...), Eqs. ( 5.16] ) remain an excellent 
approximation. [The singularity of the ZS' graph be- 
comes crucial if one is precisely interested in the value 
of r«(0 - 9') for 9 - 6' -> OB] This conclusion has been 
checked by explicit calculation at one-loop order by Chi- 
tov and SenechalE 

At each step of the RG transformation, the field is 
rescaled according to ft — [z(di)] _1 / 2 -0, where z(dt) is 
the contribution to the wave-function renormalization 
factor z(t) when the flow parameter t increases by dt: 
z(t + dt) = z(t)z(dt). This ensures that the propagator 
keeps the form (iu; — vp^k)^ 1 S3 Thus the field {jj refers 
to QP's as anticipated. 

Now we consider the renormalization of the source 
field. Although the action is expressed only in terms 
of the QP field, h (as any other external field) couples a 
priori to the bare fermions. The coupling between the ex- 
ternal field and the incoherent part of the single-particle 
spectral function appears indirectly through some renor- 
malization of the external field (see for instance Ref. |25| 
for a discussion). Here we want to eliminate such renor- 
malizations since h couples directly to the QP's. To dis- 
tinguish between the coherent (i.e., due to the QP's) and 
incoherent responses to the field h, we also consider the 



response to the field hP that couples to the incoherent 
part of the density: 



Vp 



x lim 

u 



fcf+Ao 



dk 



im U* a {k + q)^ a {k) 1. (5.29) 

:^0L ' q=Q\ 



When T 0, h n cannot create coherent particle-hole 
pairs (i.e., quasi- particlej-quasi-hole pairs) since the f2- 
limit is taken in ( |5.29 ).lij On the contrary, h couples to 
both the coherent and incoherent parts, n-loop contribu- 
tions (n > 2) to the field renormalization factor z^ h ' do 
not distinguish between h and hP . The reason is that the 
singular ZS loop does not appear in n-loop (n > 2) di- 
agrams. Therefore, the latter correspond to coupling of 
the external field to the incoherent part. On the con- 
trary, the one-loop graph considered in the preceding 
section vanishes in the f2-limit and corresponds to cou- 
pling of the external field to the coherent part. Since 
we want h to couple only to the QP's, only the one- loop 
diagram for z^ h ' has to be taken into account. QP's are 
obtained not only by filtering out the incoherent part of 
the propagator, but also by rescaling the field accord- 
ing to — [z(c?i)] _1 / 2 ^. The latter implies a concomi- 
tant rescaling of the field hi — hz(dt) to ensure that 
h^fty = h'4> * l tp' . We conclude that the only renormal- 
ization of the external field comes from the one-loop con- 
tribution to z^ considered in the preceding section. 

We also note that the diagram shown in Fig. 2b of 
Ref. [l(] is the only one contributing to the susceptibility 
in the Kadanoff- Wilson scheme, since it is the only dia- 
gram of order 0(h 2 ) (and 0(dt)) generated by the RG 
procedure. Thus, we obtain the same RG equations as in 
the one-loop calculation, the bare Fermi velocity being 
replaced by its fixed-point value v* F . Eq. (5.23) holds at 
all order in a loop expansion (when the small-angle singu- 
larity in the ZS' channel is neglected). Since the V>'s have 
been rescaled at each step of the RG transformation, we 
eventually come to 



*2-pQ* 



(5.30) 



where T^* now refers to the bare fermions, and z* is 
the fixed-point value of z(t). Eq. ( |5.30 ) agrees with the 
conclusion of Ref. 10. 



VI. CONCLUSION 

We have proposed a new microscopic description of 
Fermi liquids, which extends some early ideas of the sta- 
tistical FLT. It is based on the introduction of an effec- 
tive potential (in the sense of field theory) 5>[n], which is 
obtained from the free energy by a Legendre transforma- 
tion. 

In the more general case, the effective potential $[n] 
is a functional of the Wigncr distribution function n = 
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{^k(r(r,t)}. Small variations <$$[<5n] around the equilib- 
rium value are parameterized by the Landau / function, 
which describes the interaction between QP's. / has a 
precise microscopic definition in terms of the susceptibil- 
ity x introduced in section Using the standard as- 
sumptions of FLT, we have shown that this microscopic 
definition yields the usual identification between / and 
the f2-limit of the forward-scattering two-particle vertex. 
The effective potential 5$[<5n] yields both the response 
functions of the Fermi liquid and the quantum Boltz- 
mann equation satisfied by n\^ a (r,t). In the static and 
uniform limit, <5$[(5n] is nothing but the variation of the 
thermodynamic potential corresponding to a change Sn 
of the QP distribution function. 54>[(5n] was first intro- 
duced by Landau on phcnomcnological grounds to de- 
scribe Fermi liquids. Thus, the effective potential de- 
scribes both static and dynamic properties of Fermi liq- 
uids in a unified framework. It should be noted that this 
description does not rely on any semiclassical assump- 
tion. 

The explicit calculation of 55>[£n] can be in principle 
extended to more complicated and/or realistic situations, 
for instance by taking into account the presence of impu- 
rities and the effect of an electric field. 
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